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This  report  documents  the  results  of  a  study  aimed  at  obtaining  the 
spatial  and  temporal  distribution  of  loads  on  the  headform  as  a  result  of  a 
crown  impact  to  Army  aviator  helmets  using  a  classic  headform-helmet  drop  test. 
In  Che  actual  drop  test,  only  2  measurements  were  made:  the  impact  force  time 
history  between  the  helmet  shell  and  a  rigid  surface  and  the  acceleration  time 
history  of  the  center  of  gravity  (C.G.)  of  the  magnesium  headform. 

A  mathematical  model  of  the  drop  test  was  developed  using  a  flexible  body 
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extension  to  DADS,  a  general-purpose  nonlinear  dynastic  analysis  program 
previously  completed  at  the  University  of  Iowa,  To  exercise  the  model,  it  was 
necessary  to  experimentally  determine  the  following  helmet  and  headform 
input  parameters: 

(1)  The  center  of  gravity  (C.G.)  was  obtained  using  a  variant  of  the 
differential  weighing  technique, 

(2)  The  principal  mass  moments  of  inertia  about  the  C.G.  was  found  by  the 
repeated  application  of  the  well-known  torsional  pendulum  method. 

(3)  The  3-dimensional  (3-D)  coordinates  of  many  points  of  the  helmet  shell 
surface  was  determined  via  an  innovative  technique  for  3-D  digitization  using 
an  electronic  caliper  with  a  computer  algorithm. 

(4)  The  dimensions  of  the  gap  between  the  headform  and  helmet  liner  were 
obtained  using  an  alginate  dental  impression  material  as  a  cast  between  the 
two  objects. 

(5)  The  material  properties  of  the  helmet  shell  and  liner  were  found  through 
standard  materials  testing  methods  using  a  MTS  machine. 

In  the  mathematical  model,  the  headform  and  the  contact  surface  are 
assumed  to  be  rigid.  The  liner  material  is  modeled  as  a  series  of  nonlinear 
springs  attached  to  a  flexible  shell.  The  area  of  contact  between  the  crown 
of  the  helmet  and  the  rigid  surface  is  also  modeled  by  a  series  of  5  non¬ 
linear  springs  of  high  stiffness. 

The  model,  with  the  appropriate  geometric,  material  and  inertial 
properties  simulating  the  SPH-4  and  A1  Army  aviator  helmets,  was  run  on  the 
PRIME  750  computer.  The  results  permit  the  following  conclusions  and 
recommendations : 

I.  Comparing  the  predicted  response  of  the  SPH-4  and  the  Al,  we  see  that 
helmets  which  exhibit  similar  gross  behavior,  i.e.,  the  acceleration  of  the 
C.G.  of  the  headform  and  the  ground  contact  force,  may  have  entirely  different 
headform  contact  pressure  responses, 

II.  The  flexibility  of  the  Al  shell  allows  a  large  crown  contact  pressure 
to  develop.  The  stiffer  SPH-4  distributes  the  contact  pressure  to  the  front 
and  back  of  the  headform.  The  peak  headform  contact  pressure  is  higher  in 
the  Al  than  the  SPH-4. 

III.  For  the  liner  material  used,  a  shell  stiffness  somewhere  between  the 
Al  and  the  SPH-4  will  allow  the  most  uniform  contact  pressure  distribution. 

IV.  The  spatial  and  temporal  distribution  of  loads  on  the  rigid  headform 
could  serve  as  a  more  sensitive  method  of  comparing  the  performance  of 
different  helmets  via  the  drop  test. 

V.  The  analysis  procedure  used  here  is  capable  of  incorporating  the 
flexibility  of  the  headform,  e.g.,  humanoid  headforms  and/or  cadaveric  heads. 
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I.  Introduction 


Simulation  of  the  head  and  neck  during  an  impact  event  requires  knowledge 
of  the  spatial  and  temporal  distribution  of  loads  on  the  head.  A  headform- 
heloet  drop  test  is  often  used  to  determine  the  helmet  impact  force  and 
headform  acceleration  history.  These  global  measures,  while  useful,  are  not 
sufficient  for  a  head  and  neck  simulation.  This  study  developed  a  model  for 
helmets  that  predict  the  transient  distribution  of  loads  applied  to  the 
headform  during  an  impact  event.  In  order  to  verify  the  developed 
computational  technique,  the  headform-helmet  drop  test  is  simulated  and  the 
resulting  global  force  and  acceleration  measures  are  compared  to  experiments. 

The  Army  Aviator  production  helmet,  SPH-4 ,  aad  an  experimental  helmet 
designated  A1  were  used  in  this  study  for  both  drop  tests  and  simulations. 
Simulations  were  also  performed  for  two  helmets  with  the  same  shell  shape  and 
liner  material  as  the  SPH-4  and  A1  but  with  a  "rigid"  as  well  as  a  soft  shell. 
II.  Analysis  Procedure 

The  system  was  divided  into  four  components:  the  impact  surface,  the 
shell,  the  liner,  and  the  headform.  The  headform  and  the  impact  surface  were 
modeled  as  rigid  bodies.  Only  the  shape  and  rigid  body  inertial  properties  of 
the  headform  were  used  in  the  analysis.  Figure  1  depicts  a  plane  view  of  the 
three-dimensional  model. 

The  dynamic  analysis  was  carried  out  by  a  flexible  body  extension  to  the 
DADS^  nonlinear  dynamic  analysis  program.  This  method  for  transient  dynamic 
analysis  of  mechanical  systems  consisting  of  constrained  rigid  and  flexible 
bodies  was  developed  at  the  University  of  lowa^).  Overall  displacement  of 
elastic  bodies  was  represented  by  superposition  of  small  elastic  displacements 
on  the  large  displacement  of  the  body  reference  frames.  For  each  elastic  body 
two  seta  of  generalized  coordinates  were  employed.  First,  reference 
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generalised  coordinates  defined  the  location  of  a  body-fixed  coordinate  system 
with  respect  to  an  inertial  reference  frame*  This  coordinate  system  had  its 
origin  rigidly  affixed  to  a  point  on  the  body.  Second,  the  generalised 
coordinates  corresponding  to  elastic  deformation  were  introduced  through  a 
nodal  analysis  method.  This  set  of  coordinates  was  defined  with  respect  to 
the  body-fixed  coordinate  system. 

The  systea  equations  of  notion  were  formulated  using  Lagrange's 
equations.  Since  elastic  coordinates  associated  with  a  given  elastic  body  are 
defined  locally,  a  finite  element  approach  was  used  to  assemble  the  elements 
within  the  body,  thus  eliminating  excess  dependent  variables.  However, 
constraints  between  adjacent  bodies  were  not  used  to  explicitly  eliminate 
additional  dependent  variables.  Rather,  a  Lagrange  multiplier  technique  was 
employed.  Substructuring  methods  and  generalized  coordinate  partitioning  were 
employed  to  reduce  the  number  of  degrees  of  freedom.  The  large  displacements 
resulted  in  inertia-variant  problems  with  a  corresponding  variant  eigen- 
spectrum.  The  method  allowed  for  adjustment  of  the  number  of  elastic  degrees 
of  freedom  so  that  the  dominant  modes  are  included  in  the  dynamic  model.  The 
elastic  deformation  was  represented  by  lower  (dominant)  mode  dynamic  response 
and  higher  mode  static  effects. 

The  equations  of  motion  and  constraints  were  solved  numerically  using  a 
direct  integration  method.  Numerical  examples  have  shown  that  the  results 
obtained  by  this  method  and  those  available  in  the  literature  are  in  good 
agreement. 

As  the  deformations  of  the  helmet  shell  were  expected  to  be  small  during 
the  impact  under  consideration,  the  shell  was  modeled  using  linear  elastic 
small  deformation  plate  theory.  The  shell  is  discretized  into  three-node 
triangular  plate  finite  elements.  Figures  2  and  3  depict  the  discretized 
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model  of  the  helmet  shell.  The  SPAR  finite  element  program  was  used  to  obtain 
the  shell  stiffness  matrix,  mass  matrix,  and  mode  shapes. 

tfodeling  of  the  foam  liner  material  was  the  most  complex  aspect  of  the 
analysis.  The  material  was  expected  to  undergo  relatively  large  deformations 
during  impact  and  thus  the  response  of  the  material  cannot  be  reasonably 
modeled  by  using  a  linear  constitutive  relationship.  Since  not  all  of  the 
liner  is  in  initial  contact  with  the  headform,  the  gap  between  the  headform 
and  the  liner  material  must  be  taken  into  account  in  the  simulation.  In  order 
to  model  these  aspects  of  the  liner,  it  was  divided  into  small  regions  and 
each  region  is  represented  by  a  nonlinear  one-dimensional  element.  A  typical 
force-deformation  curve  for  an  element  is  shown  in  Fig.  4.  The  gan  was  taken 
to  be  the  average  distance  between  the  liner  and  the  headform  over  the  element 
area.  Until  the  headform  has  advanced  to  close  the  gap,  there  was  no  response 
from  the  liner.  The  use  of  a  trilinear  force-deformation  relation  was  well 
justified  by  the  form  of  the  experimentally  determined  force-deformation 
relation  shown  in  Fig.  5.  Several  typical  liner  materials  have  been  observed 
to  display  this  stiffening  character. 

The  impact  surface  was  modeled  as  a  set  of  linear  elements  with  a  gap 
capability.  These  elements  were  distributed  around  the  crown  of  the  helmet 
over  the  region  which  was  reasonably  expected  to  come  into  contact  with  the 
impact  surface.  The  impact  surface  was  assumed  to  be  flat  and  the  element 
gaps  were  determined  by  the  shell  geometry. 

The  initial  condition  assumed  the  same  pre-impact  velocity  to  the  helmet 
and  headform.  The  headform  initial  position  was  one  that  just  touched  the 
helmet  liner  at  two  points  located  in  the  posterior  and  anterior  regions  of 
the  plane  of  mid-sagittal  symmetry.  An  initial  gap  between  the  helmet  liner 
and  headform  existed  over  the  crown  region.  The  initial  position  of  the 
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helmet  was  one  that  yielded  a  small  initial  gap  between  the  shell  and  the 
impact  surface,  i.e.,  the  condition  of  the  system  at  the  instant  before  the 
impact. 

A  constant  gravitational  force  was  included  on  the  helmet  and  the 
headform.  The  helmet  was  constrained  to  move  straight  downward  and  rotation 
was  prevented.  The  headform  was  also  constrained  to  move  straight  downward 
but  rotations  were  allowed. 

In  this  study  the  simulation  was  performed  from  the  instant  before  impact 
until  the  time  when  the  helmet  has  rebounded  from  the  impact  surface. 

III.  Material  and  Geometric  Properties 

A.  Shell 

The  shell  material  properties  for  the  A1  and  SPH-4  helmets  were  obtained 
by  material  testing.  All  samples  were  cut  with  a  table  saw  using  a  laminated 
blade.  Hie  width  of  the  specimens  was  approximately  1  cm.  for  a  length  of  8 
to  10  cm.  The  samples  from  the  Army  helmets  were  cut  such  that  half  of  them 
had  a  surface  layer  fiber  orientation  directed  along  the  length  of  the 
specimen.  The  remainder  were  cut  such  that  the  surface  fiber  direction  was 
transverse  to  the  axis  of  loading.  A  total  of  15  tests  were  conducted,  three 
for  each  condition.  For  typical  results,  see  Figures  6  and  7.  Density  was 
determined  using  a  volumetric  displacement  procedure. 

All  tests  were  conducted  on  a  closed  loop  servo-hydraulic  Materials 
Testing  System  (MTS)  machine.  Specimens  were  mounted  in  friction  grips  and 
tests  were  run  with  a  ramp  extension  rate  of  1.27  mm/min  (0.05  in/min). 

Strain  was  measured  using  120  fl  strain  gages  in  a  Wheatstone  bridge  circuit. 
Two  arm  bridges  were  used  for  both  longitudinal  and  latitudinal  strains  and  so 
placed  that  bending  effects  were  eliminated.  The  bridge  circuit  was  amplified 
using  Honeywell  Accudata  218  signal  conditioners.  Hardcopy  output  was 
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recorded  on  an  X-Y  plotter  which  graphed  the  force  output  versus  the  strain, 
obtained • 

Due  to  the  strip  nature  of  the  specimens,  no  complete  stress-strain 
curves  were  generated.  However,  for  the  purposes  of  this  study.  Young's 
modulus,  E,  and  Poisson's  ratio,  u  ,  were  the  only  material  properties 
needed.  Some  difficulty  was  encountered  while  testing  the  Army  A1  helmet 
specimens.  At  approximately  250  N  the  fiberglass  layers  failed  in  shear  and 
the  specimens  would  pull  out  of  the  grips.  Until  this  point  though,  there  was 
little  or  no  relative  displacement  between  the  fiberglass  layers.  This  was 
determined  by  coating  the  edges  with  a  polyurethane  layer,  step  loading  the 
specimen,  and  at  each  load  increment,  looking  for  cracks  in  the  plastic.  As 
can  be  seen  in  Fig.  6  there  was  little  difference  in  the  curves  with  respect 
to  different  surface  fiber  directions.  For  the  modeling  aspect  of  this  report 
a  single  Young's  modulus  was  chosen  for  each  Army  helmet. 

B.  Foam  Liner 

New  helmet  foam  liners  and  the  liners  from  an  A1  and  SPH-4  were  tested. 
The  density  of  all  liners  was  determined  using  a  volumetric  displacement 
procedure  and  was  found  to  be  0.077  *  0.003  g/cm.  No  difference  in  impact 
characteristics  were  3een  between  different  liner  sources  and  they  will  not  be 
distinguished  in  the  remainder  of  this  report.  Sheet  specimens  were  taken 
from  the  lateral  aspect  of  the  crown  such  that  che  radius  of  curvature  was 
minimized.  Cylindrical  samples  were  chosen  at  random  and  were  cut  with  a  one 
inch  plug  cutter. 

All  impact  tests  were  conducted  on  a  closed  loop  servo-hydraulic 
Materials  Testing  System  (MTS)  machine  (Model  810)  rated  at  10  kN. 
Displacements  were  monitored  with  a  LVDT  attached  as  an  integral  part  of  the 
MTC  actuator.  Hardcopy  output  was  a  photograph  from  a  storage  oscilloscope 
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which  displayed  the  outputs  of  the  LVDT  and  MIS  machine's  load  cell. 
Displacement  control,  with  a  triangle  wave  forcing  function,  was  used  and  peak 
displacements  ranged  from  3  to  12  am.  Rise  times  varied  from  12  to  26  msec. 
All  tests  were  run  with  the  specimen  being  placed  on  a  large  steel  plate 
bolted  to  the  top  of  the  actuator.  Specimens  were  oriented  such  that  the  line 
of  impact  would  be  normal  to  the  inner  liner  surface.  The  impactor  was  a  one- 
inch  diameter  steel  rod  attached  to  a  stationary  load  cell  (see  Fig.  8). 
Initial  geometries  were  determined  with  the  aid  of  a  precision  micrometer. 

As  indicated  above,  two  different  types  of  tests  were  run.  The  first  of 
these  was  with  a  piece  of  foam  liner  where  the  area  was  much  larger  than  that 
of  the  impactor,  i.e.,  an  “infinite"  sheet.  Due  to  the  radius  of  curvature  a 
backing  was  needed  for  complete  support.  This  was  achieved  by  molding  PMMA  to 
the  piece  of  foam  in  the  orientation  desired.  A  preload  of  5  N  was  used  to 
maintain  position  prior  to  the  test.  The  second  type  of  specimen  was 
cylindrical  in  shape  with  the  same  diameter  as  that  of  the  impactor.  These 
tests  were  run  with  a  20  N  preload.  As  would  be  expected,  the  two  different 
types  of  tests  had  similar  shaped  force  displacoent  curves.  However,  sheet 
specimens  tended  to  absorb  approximately  15  percent  more  energy  due  to  tearing 
of  the  compressed  plug  from  the  surrounding  material  (see  Figures  9-11).  All 
curves  were  characterized  by  an  initial  linear  region  where  the  foam 
approximated  an  isotropic  linear  elastic  material  (see  curve  7  in  Figure 
11).  At  about  20  percent  strain  and  500  N  there  was  a  fairly  abrupt  region  of 
transition.  Beyond  this  point  there  was  a  much  lower  slope  where  disruption 
of  the  foam  cellular  structure  occurs.  It  was  in  this  area  that  most  of  the 
energy  absorption  took  plac-».  There  was  a  second  transition  region  (see  curve 
3  in  Fig.  9)  at  about  75  percent  strain  where  the  stiffness  increased.  This 
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behavior  has  been  termed  "lock-up"  and  occurred  when  the  cellular  form  of  the 
foaa  acted  as  a  stiff  unit*  At  this  strain  level  the  foam  was  almost  a 
uniformly  dense  material  with  very  little  energy  absorption  capabilities.  It 
should  be  kept  in  stirtd  that  the  impact  force-displacement  curves  were  very 
much  dependent  on  loading  races,  i.e.,  rise  times. 

C.  Gap  Generation 

The  contact  area  between  the  metal  headfora  and  the  helmet  liner  was 
experimentally  determined.  Selected  points  on  the  helmet  shell  were  marked. 

At  these  points  1/16  inch  holes  were  drilled  through  the  helmet  shell  and  the 
liner.  This  hole  was  angled  such  that  it  would  be  normal  to  the  shell 
surface.  Using  a  pen,  these  points  were  labeled  by  indenting  the  foam 
liner.  The  metal  headform  was  also  marked  using  labeling  tape  (periods  for 
the  points)  with  an  identifying  numerical  tag. 

An  alginate  dental  impression  material  Jeltrate®  was  used  for  casting  the 
gap.  This  material  when  mixed  is  a  liquid  of  low  viscosity  but  quickly  (180 
sec)  sets  to  a  rubbery  consistency.  The  helmet  was  filled  with  the  liquid 
impression  material  and  the  metal  headform  positioned.  Once  the  Jeltrate®  had 
set,  the  headform  and  cast  were  removed.  The  casting  clearly  showed  the  metal 
headform  and  the  helmet  liner  points.  The  cast  was  then  cut  at  these  points 
and  the  thickness  measured  with  a  micrometer. 

For  several  selected  helmet  drop  tests  a  sheet  of  carbon  paper  was  placed 
between  the  metal  headform  and  the  helmet  foam  liner.  The  carbon  side  was 
positioned  facing  the  liner.  The  drop  tests  were  then  conducted  according  to 
standard  procedure.  In  this  manner  regions  vere  marked  on  the  liner  where  the 
headform  caused  deformations  during  the  drop.  In  addition,  using  shading 
patterns,  rough  estimates  of  relative  displacements  could  be  made.  The 
outlines  of  these  impact  contact  regions  were  demarcated  by  indenting  the 
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liner  with  a  pen.  Then  using  the  same  procedure  as  above  a  casting  was 
made.  The  cast  was  highlighted  with  ink  and  then  photographed. 

IV.  Helmet  —  Bead form  Drop  Tests 

The  behcvior  of  the  helmet  in  a  dynamic  environment  was  evaluated  by  drop 
tests.  A  medium-size  magnesium  headform  was  provided  by  the  Contract 
technical  Monitor  (CTM),  U.S.  Army  Aeromedical  Research  Laboratories,  Ft. 
Rucker,  AL.  At  the  approximate  center  of  gravity  of  the  headform  an 
accelerometer  (Endevco  2264,  2000  G)  was  attached  with  its  sensitive  axis 
oriented  vertically-  The  drop  method  conformed  in  general  to  the  ANSI/ASTM  F 
429-75  standard.  The  testing  setup,  purchased  from  Humanoid  Systems,  Canson, 
CA,  used  no  guide  assembly  since  only  crown  impacts  were  conducted.  The 
experimental,  setup  shown  in  Fig.  12  consisted  of  a  channel  frame  with  an 
electromagnet  attached  at  the  center  of  its  horizontal  cross  member.  The 
helmet,  fitted  with  the  headform,  was  held  by  the  electromagnet.  The  base  of 
the  apparatus  was  a  structural  steel  plate  5.08  cm  thick.  The  midportion  of 
the  plate  was  machined  out  to  accomodate  CNSTUL  20  kN  force  transducer 
(manufactured  by  PCB  Piezotronics,  Inc.).  The  height  from  the  base  to  the 
electromagnet  was  86  cm.  For  both  the  SPH-4  and  the  A1  helmet,  the  drop 
height  was  approximately  59  cm. 

The  drop  test  begins  with  the  interruption  of  the  current  to  the 
electromagnet  by  flipping  a  switch.  The  transient  events  were  recorded  by  the 
appropriate  above-mentioned  transducers,  signal  conditioners  (Honeywell 
Accudata  218  series),  tape  recorder  (Honeywell  FM  Tape  Recorder  System  5600E 
and  storage  oscilloscope  (Tektronix  7623A).  The  acquired  analog  signals  were 
converted  to  digital  form,  through  A/D  conversion  available  on  our  PDP  11/34 
minicomputer  for  further  processing.  The  final  display  can  be  either  in 
Polaroid  photographs  of  the  storage  oscilloscope  or  computer  graphics  display. 
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V,  Models 


Initially  three  separate  simulations  were  performed:  one  with  a  rigid 
helaet  shell,  one  with  a  helmet  shell  model  based  on  the  Al,  and  one  with  a 
model  based  on  the  SPH-4.  The  shapes  of  the  helmet  shell  and  liner  were  found 
to  be  identical  in  the  Al  and  the  SPH-4.  The  shell  material  in  the  A1  is  much 

less  stiff  than  that  of  the  SPH-4  and  the  Al  material  is  thinner.  The  three 

simulations  represented  a  series  of  helmets  with  identical  shell  shape  and 
liner  material  but  differing  shell  stiffness. 

The  shell  was  constructed  of  a  laminated  fibrous  composite  material.  The 
orientations  of  the  lamina  varied  from  point  to  point  on  the  shell  and  from 
helmet  to  helmet  at  the  same  point.  A  number  of  rectangular  specimens  were 
cut  from  each  helmet  at  various  locations  and  orientations.  It  was  found  that 
the  material  properties  of  the  laminate  are  not  heavily  orientation 
dependent.  The  material  was  then  assumed  to  be  transversely  isotropic.  The 

values  obtained  in  the  tensile  tests  were  averaged  to  obtain  a  single  Young’s 

modulus  and  Poissoa  ratio  for  each  shell  material.  Shell  material  properties 
used  in  the  analysis  are  presented  in  Table  l. 

Sample  helmets  were  cut  into  radial  sections  so  that  the  shell  and  liner 
thickness  could  be  measured.  Neither  the  shell  nor  the  liner  thickness  varied 
abruptly.  The  shell  shape  and  the  headform  shape  were  measured  using  either  a 
3-D  sonic  digitizer  (Graf-Pen)  or  the  electronic  caliper  system  described  in 
Appendix  A.  The  measured  shell  shape  was  used  to  construct  the  finite  element 
model  mesh. 

It  was  expected  that  the  contact  stress  distribution  will  depend  heavily 
on  the  initial  gap  between  the  headform  and  liner.  Rather  than  computing  this 
gap  by  differencing  the  measured  headform  and  shell  plus  liner  shape,  it  was 
decided  to  measure  the  gap  directly.  A  quick  setting  liquid  plastic  was 


poured  into  an  Inverted  helmet  with  the  headform  In  place.  After  the  plastic 
had  set,  it  was  removed,  sectioned  and  Its  thickness  measured.  These 
measurements  were  used  to  specify  the  initial  gap  used  in  the  analysis. 

The  S?AB.(3)  finite  element  program  was  used  to  compute  the  shell 
stiffness  matrix,  mass  matrix,  eigenvalues,  and  mode  shapes.  The  mass  of  the 
liner  material  was  distributed  over  the  shell  in  the  areas  of  labeled  nodes  in 
Fig.  13.  Each  finite  element  shell  model  consisted  of  83  nodes  connected  by 
triangular  plate  elements.  The  first  16  modes  were  used  In  the  analysis  of 
the  SPH-4  and  the  first  17  modes  were  used  for  the  Al.  Table  2  contains  the 
eigenvalues  used.  Plots  of  the  mode  shapes  are  found  in  Appendix  B. 

Compression  tests  were  performed  on  samples  of  the  shell  liner  material. 
Figure  5  displays  some  typical  results.  Based  on  these  tests  the  liner 
material  was  modeled  as  being  trilinear  as  shown  in  Fig.  4.  The  stress-strain 
slopes  used  are  found  in  Table  3.  The  effective  spring  constants  of  the  one- 
dioensional  elements  used  to  model  the  liner  were  computed  on  the  basis  of  the 
material  properties  and  the  extent  of  liner  material  around  each  node.  A 
viscous  damping  factor  was  included  in  the  liner  material. 

The  liner  was  modeled  as  32  one-dimensional  elements  that  were  attached 
to  the  shell  at  the  points  labeled  in  Fig.  13.  The  area  over  which  the  liner 
elements  are  distributed  was  chosen  to  cover  the  possible  contact  3rea  noted 
from  the  carbon  paper  impressions  made  during  drop  tests.  An  additional  5 
linear  one-dimensional  gap  elements  were  used  tc  model  the  impact  surface. 
These  were  attached  to  the  nodes  labeled  3,  12,  13,  29  and  64  in  Fig.  13. 

VI.  Results  and  Discussion 

The  headform  helmet  drop  test  was  simulated  using  the  finite  element 
model  of  Army  helmets  developed  in  Section  IV.  In  this  section,  computational 
results  for  the  SPH-4,  Al  and  the  rigid  shell  are  presented. 
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Figures  14(a)  to  (d)  are  the  plots  of  the  numerical  predictions  for  the 
SPH-4  helmet,  corresponding  to  (a)  position,  (b)  velocity,  (c)  acceleration, 
and  (d)  impact  force  between  the  ground  and  the  SPH-4  shell.  Similarly,  Figs. 
15(a-d)  and  I6(a-d)  present  the  corresponding  results  for  the  Al  and  rigid 
shell,  respectively.  In  Fig.  14(a),  the  solid  line  shows  the  position  of  the 
crown  of  the  SPH-4  shell  (node  3  in  Fig.  3)  and  the  dotted  line  shows  the 
position  of  the  crown  of  the  headform,  as  functions  of  time.  Nonlinear 
springs,  simulating  the  ground,  became  active  when  the  distance  between  the 
origin  and  the  shell  surface  was  less  than  10  millimeters.  In  these  figures, 
it  was  observed  that  the  distance  between  the  SPH-4  shell  and  the  headform  did 
not  change  very  much  during  impact.  The  same  was  true  for  the  rigid  shell,  as 
shown  in  Fig.  16(a).  On  the  other  hand,  the  distance  between  the  Al  shell  and 
the  headform  changes  substantially  as  shown  in  Fig.  15(a).  This  was  an 
expected  result,  since  the  stiffness  of  the  Al  shell  was  about  30  times  lower 
than  that  for  the  SPH-4.  The  Al  shell  can  deform  more  easily  so  that  two 
points  can  approach  more  closely  than  those  in  the  SPH-4.  Figure  14(b)  shows 
the  velocity  history  of  the  SPH-4  shell  and  the  headform.  During  experimental 
head  drop,  the  velocity  of  the  shell  and  headform  at  impact  was  -3.42  o/sec. 
Therefore,  the  initial  conditions  of  the  numerical  experiment  were  chosen  so 
that  the  same  velocity  can  be  observed.  In  Fig.  14(b),  both  the  SPH-4  shell 
and  the  headform  approached  the  ground  at  a  downward  velocity  of  -3.24 
m/sec.  The  shell,  hitting  the  ground  first  obtained  positive  upward  velocity 
gradually  after  t  ■  0.004  sec.  due  to  the  upward  reaction  force  from  the 
ground.  The  headform  moving  toward  the  ground,  hits  the  liner  at  t  *  0.006 
sec  giving  the  shell  a  negative  velocity.  The  shell  bounces  between  the 
ground  surface  and  the  headform.  This  oscillatory  behavior  can  be  observed  in 
Fig.  14(b).  The  numerical  solution  for  the  rigid  shell  gave  almost  identical 


results  as  that  of  the  SPH-4,  as  shown  in  Fig.  16(b).  Because  of  the 
flexibility  of  the  shell,  the  numerical  results  for  the  A1  helmet  gave  a 
smaller  amplitude  of  velocity  oscillatory  behavior  (Fig.  15(b})  than  the  ocher 
two  helmets. 

Fig.  14(c)  is  a  plot  of  the  acceleration  history  of  the  SPH-4  shell  and 
the  headform.  In  this  figure,  one  can  see  the  reaction  forces  between  the 
shell  and  the  headform  as  a  result  of  impact.  At  t  ■  0.008  sec.,  the  headform 
obtains  an  upward  acceleration  and  the  shell  receives  downward  acceleration 
due  to  the  impact  between  these  two.  By  observing  this  feature,  one  can  see 
that  the  next  impact  happens  just  after  0.01  sec.  Thus,  the  acceleration  of 
the  shell  changes  sign  several  times  during  a  single  impact  event.  The 
acceleration  history  of  the  rigid  shell  yields  a  similar  magnitude  and  pattern 
but  fewer  impacts  between  the  shell  and  the  headform  can  be  seen  in  Fig. 

16(c).  In  Fig.  14(c),  the  acceleration  of  the  headform  for  the  SPH-4  shell 
has  two  peaks.  One  peak  is  1900  (m/sec^)  *  194  g  at  t  *  0.008  and  the  second 
is  1250  (m/sec^)  >  127  g  at  t  -  0.015,  with  a  time  difference  of  0.0025  sec. 
Figure  17  is  a  plot  of  the  experimental  results.  The  observed  acceleration 
history  of  the  headform  also  has  two  peaks.  The  first  peak  has  a  magnitude  of 
160  g  and  the  second  peak  130  g,  with  a  time  difference  of  0.002  sec.  Thus, 
the  acceleration  histories  of  the  numerical  simulation  and  the  actual 
experiment  are  in  good  agreement,  not  only  in  pattern  but  also  in  their 
magnitudes. 

Figure  15(c)  shows  the  numerical  result  of  the  acceleration  time 
histories  for  the  A1  headform  and  helmet.  One  observes  that  the  peak  headform 
acceleration  has  a  magnitude  1300  (m/sec^)  ■  132  g.  Figure  18  is  a  plat  of 
the  acceleration  of  the  Al  experimental  data,  whose  peak  is  120  g.  Though  the 
patterns  look  different,  the  magnitudes  are  close. 
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Figure  14(d)  shows  the  computed  impact  force  between  the  SPH-4  shell  and 
the  ground,  and  Fig.  17  shows  the  corresponding  experimental  result.  For  the 
computed  results,  the  peak  force  is  3700  N  with  a  time  duration  of  0.011  sec., 
while  the  experimental  results  show  a  peak  force  of  7500  N  with  a  time 
duration  of  0.010  sec.  Even  though  the  computed  force  is  almost  half  the 
magnitude  of  the  experimental  results,  the  patterns  are  similar,  having  two 
almost  Identical  peaks  and  time  durations.  Computed  results  for  the  rigid 
shell  are  plotted  in  Fig.  16(d)  showing  a  peak  value  slightly  higher  than  the 
results  for  the  SPH-4.  Figure  15(d)  is  a  plot  of  the  computed  impact  force 
between  the  A1  helmet  and  the  ground  and  Fig.  18  shows  the  experimental 
results  for  the  Al  helmet.  The  peak  force  of  the  computed  results  was  3700  N 
and  for  the  experimental  results,  6000  N.  The  pattern  is  also  different  and 
the  numerical  results  has  two  peaks  while  the  experimental  results  are 
unimodal. 

Figures  19  (a)  to  (c)  are  the  distribution  history  for  the  pressure 
between  the  shell  and  the  headform  for  the  SPB-4  helmet.  First,  the  spatial 
three-dimensional  helmet  is  projected  to  the  x-y  plane.  Then  the  pressure 
levels  at  each  node  are  expressed  as  the  z-coordinate.  Only  selected  time 
instances  are  presented  here.  At  each  time  instance,  two  figures  are  drawn. 
One  is  a  side  view  and  the  other  a  rear  view  of  the  pressure  distribution.  At 
t  ■  0.0039  sec.,  no  pressure  is  observed  between  the  shell  and  the  headform. 

At  t  *  0.0051  sec.,  some  small  pressures  are  observed  at  the  front  and  rear 
portions  of  the  helmet.  At  time  t  *  0.0060  sec.,  shown  in  Fig.  19(a),  no 
significant  contact  pressure  is  observed,  even  though  the  shell  is  at  the 
lowest  position  in  Fig.  14(a),  the  impact  between  shell  and  headform  has  not 
occurred  yet.  At  time  t  *  0.0072,  some  of  the  liner  material  is  compressed 
due  to  the  impact  and  high  pressure  observed  at  the  rear  portion  of  the  shell 
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(Node  10)  with  values  of  0.58  x  10^  N/®"”  at  time  t  *  0.0072.  Similar  high 
pressure  appears  at  the  front  portion  of  the  shell  (node  4)  at  time  t  - 
0.0075.  The  contact  area  increases  at  time  t  *  G.0081.  Then  after  t  * 

0.0099,  those  contact  pressures  decrease  and  die  out  at  t  -  0.0159  sec. 

During  the  impact,  the  crown  of  the  shell  (Node  3)  never  experiences  the 
contact  pressure.  This  is  in  good  agreement  with  the  experimental  result. 

Figures  20(a)  to  (c)  show  the  pressure  distribution  of  the  Al  helmet. 

The  behavior  is  very  different  from  the  SPH-4  and  it  has  some  pressure  at  the 
crown  (Node  3)  at  time  t  *  0.0069  and  t  -  0.0072.  At  time  t  -  0.0075  sec., 
high  pressure  appears  at  the  rear  portion  (Node  10)  and  very  high  pressure 
(0.52  x  106  N/m^)  is  generated  at  the  crown  at  time  t  ■»  0.0081  sec.  The  same 
pressure  can  be  observed  at  t  *  0.0099  sec.  After  t  ■  0.012  sec.,  the 
pressure  decreases  at  the  crown  and  the  pressures  at  the  front  and  rear 
portions  are  not  significant.  The  pressure  at  the  crown  disappears  at  t  * 
0.0159  sec.  The  significant  difference  between  SPH-4  and  Al  is  that  the  front 
and  rear  portions  are  load  bearing  portions  for  SPH-4,  while  the  crown  is  the 
load  bearing  part  for  the  Al  helmet. 

Figures  21(a)  to  (c)  show  pressure  distribution  for  the  rigid  shell  and 
gives  behavior  similar  to  the  SPH-4  helmet,  i.e.,  the  SPH-4  helmet  is  more 
like  a  rigid  shell. 

Although  the  global  SPH-4  predictions  agree  well  with  the  drop  test 
results,  the  agreement  of  the  Al  predictions  with  the  drop  test  is  less 
satisfactory.  In  order  to  further  understand  the  effect  of  shell  stiffness  on 
the  global  predictions,  a  "soft”  shell  model  was  constructed.  The  shell  of 
this  model  again  used  the  same  geometry  as  the  SPH-4  and  the  Al,  but  the  shell 
stiffness  was  30  times  smaller  than  that  of  the  Al.  The  liner  material 
characteristics  and  distribution  were  identical  to  those  used  in  the  Al  and 
SPH-4. 
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Figures  22  (a)  to  (d)  are  plots  of  the  predicted  position,  velocity, 
acceleration  and  impact  force  between  the  "soft"  shell  and  ground.  The  impact 
force  prediction  for  this  model  agrees  qualitatively  with  that  seen 
experimentally  for  the  Al  shell.  The  two  peaks  for  the  force  curve  (see  the 
SPH-4  case)  have  merged  into  a  single  peak. 


VIZ.  Conclusions 


From  the  above  results,  one  can  conclude  the  following: 

1.  The  analysis  procedure  used  is  capable  of  good  qualitative 
agreement  with  experimental  results.  Features  observed  in  the 
gross  behavior,  i.e.,  acceleration  of  the  center  of  gravity  of 
the  headforo  and  the  ground  contact  force,  are  predicted. 

2.  From  the  analysis,  it  is  confirmed  that  the  multiple  peaks  seen 
in  the  contact  force  in  drop  tests  are  due  to  the  shell  bouncing 
between  the  headform  and  the  ground  during  the  contact  event. 

3.  Comparing  the  predicted  response  of  the  SPH-4  and  the  A1  we  see 
that  helmets  which  exhibit  similar  gross  behavior  may  have 
entirely  different  headform  contact  pressure  responses. 

4.  Helmet  shell  flexibility  can  affect  gross  response  measures  only 
for  a  sufficiently  large  change  in  shell  stiffness. 

5.  Comparing  the  predicted  headform  contact  pressure  distribution 
of  the  SPH-4  and  the  A1  we  note  that  the  flexibility  of  the  A1 
allows  a  large  crown  contact  pressure  to  develop.  The  stiffer 
SPH-4  distributes  the  contact  pressure  over  two  regions  in  the 
front  and  back  of  the  headform.  The  peak  headform  contact 
pressure  is  higher  for  the  A1  than  for  the  SPH-4. 

6.  The  headform  acceleration  and  the  shell  contact  force  time 
history  are  gross  measures  of  the  helmet  performance.  Whether 
or  not  the  pressure  magnitudes  are  diminished,  and  by  how  much, 
can  only  be  ascertained  by  examining  the  details  of  the  contact 
impact  phenomenon. 


18 


VIU .  Recommendations 


l*  the  most  uniform  contact  pressure  distribution  is  produced  by  a 
shell  that  is  sufficiently  stiff  to  prevent  "large"  local 
deformations  at  the  point  of  impact,  yet  sufficiently  flexible  to 
allow  a  wide  load  distribution.  For  the  liner  material  studied,  a 
shell  stiffness  somewhere  between  that  of  the  A1  and  the  SPH-4  is 
recommended.  It  is  suggested  that  a  parametric  optimization  be 
conducted  to  determine  the  optimal  shell  stiffness. 

2.  Increasing  the  liner  thickness  in  general  decreases  the  initial 
stiffness  but  may  also  shift  the  transition  or  "toe"  regions  of 
the  load-deformation  curve.  The  effects  of  increasing  liner 
thickness  should  be  experimentally  determined  prior  to  exercising 
the  model  for  different  liners  and/or  thicknesses. 

3.  To  augment  the  headform-helmet  drop  tests,  ensemble  quasistatic 
and  dynamic  tests  of  the  helnet-headfora  system  in  a  servo- 
controlled  closed  loop  hydraulic  material  test  machine  should  be 
undertaken. 

4.  The  flexibility  of  the  headform  in  drop  tests  involving  humanoid 
headforms  and  cadavera  should  be  incorporated  in  any  future 
analyses.  The  analysis  procedure  used  here  is  potentially  capable 
of  performing  such  a  task;  however,  a  general-purpose  finite 
element  program,  e.g.,  ANSYS,  will  be  less  time  consuming  in  the 
long  run. 

5.  The  spatial  and  temporal  distribution  of  loads  on  the  head 
predicted  by  this  model  could  serve  as  a  more  appropriate  load 
input  to  the  homeomorphic  finite  element  model  of  the  head  and 
neck^) . 
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Young’s  modulus  E 
Poison's  ratio  v 
Hass  density  p 


Table  t:  Shell  Material  Properties 


0.6889x10^  erg/msec3 


Eigenvalues : 


ciwui&A  jsafc.  iiwt.  *!UE  i/rW*  fi* 


17.914x10^  erg/msec^ 
0.25 

2100.0  erg/m3 


2100.0 


0.22 

erg/m3 


Table  2:  Helmet  Eigenvalues 


1) 

6. 902x1 03 

8.912x10* 

2) 

1.118xl04 

9.056x10* 

3) 

1.942xl05 

9.l57xl6 

4) 

3.1068xl05 

9.235x10* 

5) 

6.724xl05 

9.380x10* 

6) 

1.9150x10® 

9.562x10* 

7) 

4.0013x10® 

9.607x10* 

8) 

4.9206x10® 

9.734x10' 

9) 

5.7738x10® 

9.744x10* 

10) 

6.0937x10® 

10.04x10® 

ID 

9.845x10® 

10.23x10® 

12) 

1.275xl07 

10.39xLG® 

13) 

1.8587xl07 

10.42x10® 

14) 

2.1561xl07 

10.49x10® 

15) 

2.6467xl07 

10.77x10® 

16) 

2. 6842xl07 

10. T  7xl0® 

17) 

11.21x10® 
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Table  3:  Liner  Element  Characteristics 


Spring  Displacement  <  0  0-3  3-7  7-15 


Stiffness  Modulus  0  5  x  106  1  a  106  5  x  107 

(Pa) 


Displacement 
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Fig.  14(c)  Time  History  of  the  Acceleration  of 

Shell  and  He ad form  at  the  Crowns  for  5PH4 

38 


Time  (sec) 


Fig.  15(a)  Time  History  of  the  Position  of  Crowns 
of  Shell  and  Headform  for  A1  Helmet 


tion  (ro/aec 


Time  (sec) 


Fig.  3.6(a)  Time  History  of  the  Position  of  Crowns 
of  Shell  and  Headform  for  Rigid  Helmet 


Force  (newton) 


Rigid  Helmet 


Time  ( sec ) 

Fig.  16(d)  Time  History  of  Impact  Force  between 
Shell  and  Ground  for  Rigid  Helmet 


Pig.  18  Experimental  Results  of  Force  and  Acceleration  for  A1 
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Fig*  19  (c)  Pressure  Distribution  of  SPH4 
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The  inertia!  and  geometrical  properties  of 
helmets 

GLEN  O.  NJUS,  Y.  KING  LIU.  aad  THOMAS  A.  NYE 

Department  of  Biomedical  Engineering,  College  of  Engineering,  The  University  of  Iowa, 

Iowa  City,  IA  52242 


abstract 

NJUS.  GLEN  o..  Y.  KING  UU.  tad  THOMAS  A.  NYE.  The  in- 
otia!  sad  geocoetricel  properties  of  helmets.  Med.  Set  Sports  Ex- 
ere..  Vet  16.  No.  5.  pp.  494-505.  1984.  The  center  of  prsvity 
(CSC)  aad  the  pnacipel  anas  moments  of  inert!*  tbout  the  CG  of 
Amy  tvistor,  America*  football  aad  bicycle  helmet*  were  ex¬ 
perimentally  determined  by  a  vtnetton  of  the  rlmac  dt&ereaBai 
weighing  and  torsional  peodktWi  techniques,  la  the  course  of 
these  experimmtt*.  an  innovative  method  for  threeshtnensonal 
(3D)  iftgMtsstton  was  fotndL  An  electronic  caliper,  which  mea¬ 
sured  length,  was  used  with  a  computer  algorithm  to  achieve  30 
dlgitixstkm.  The  results  of  the  above  measarmneats  show  that  the 
weight  of  the  helmet  and  the  distances  iron  die  OG  to  the  or¬ 
thogonal  coordinate  axes  intercepts  with  the  outer  shell  surface 
were  highly  correlated  with  its  principal  mass  moments  of  inertia 
A  set  of  regression  equations  was  derived  on  theoretical  consid¬ 
erations  aad  served  to  unify  the  experimentally  obtained  data. 
Our  -multi  indicate  that  the  principal  mats  moments  of  taenia  of 
helmets  vary  linearly  with  its  mass  but  noolinearly  with  dzt  and 
shape.  For  a  hairnet,  given  its  weight  and  certain  geometrical 
distances,  the  regression  equations  estimate  the  principal  mass 
ntomcrti  of  inertia  to  within  5%  of  its  experimentally  determined 
value*.  For  the  helmets  studied  in  this  aeries,  a  modified  linear- 
regresdon  relationship  between  the  principal  mass  moments  of 
inertia  and  its  mass  was  found.  This  result  is  reasonable  because 
the  mass  distribution  of  the  current  generation  of  helmets  arc  set 
primarily  by  the  head  dxe  and  secondarily  by  helmet  size,  shape. 
awl  materials. 

FOOTBALL  HELMET.  AVUTOH  HELMET.  CENTER  OF  GRAVITY, 

MASS  MOMENT  OF  INERTIA.  BIOMECHANICS 


The  determination  of  the  center  of  gravity  (CG)  and 
the  magnitude  of  the  principal  mass  moments  of  inertia 
(with  its  assviated  principal  directions)  arre  prerequi¬ 
sites  to  the  solution  of  any  problem  in  rigid-body  dynam¬ 
ics.  Several  experimental  methods  exist  to  determine  these 
inertial  properties  of  a  rigid  body. 

For  a  Sat  object,  the  method  of  suspension  is  ideal.  Hie 
object  is  placed  in  suspension  about  some  arbitrary  point 
until  equilibrium,  then  the  center  of  gravity  must  He  on 
a  vertical  line  through  the  point  of  suspension;  the  pro¬ 
cedure  is  repeated  for  a  second  arbitrary  point.  The  cen¬ 
ter  o*  gravity  lies  at  the  point  of  intersection  of  the  first 
and  second  lines.  The  method  of  balance  consists  of  rest¬ 
ing  the  fiat  object  on  a  knife  edge  until  precarious  equi¬ 
librium  is  achieved;  the  procedure  is  also  repeated  for 
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additional  orientations.  The  intersection  of  the  lines 
formed  by  the  knife  edge  determines  the  center  of  grav¬ 
ity. 

For  an  irregular  three-dimensional  (3D)  object,  the 
above  procedure  needs  considerable  modification  before 
it  becomes  practical.  A  method  developed  by  Williams 
and  Lissner  (5)  for  finding  the  weight  of  an  in  vivo  body 
segment  consists  of  weighing  the  entire  body  on  two  knife 
edges,  with  one  of  the  knife  edges  resting  on  a  weight 
scale.  The  distances  from  the  other  knife  edge  to  the 
centers  of  gravity  are  assumed  known.  The  center  of  grav¬ 
ity  of  the  segment  of  interest  is  altered  through  an  ap¬ 
propriate  joint.  The  changes  in  the  center  of  gravity  lo¬ 
cation  and  the  change  in  scale  readings  between  the  first 
and  second  weighing  are  noted-  The  segment  weight  is 
equal  to  the  distance  between  the  knife  edges  multiplied 
by  the  difference  between  scale  readings  and  divided  by 
the  horizontal  displacement  of  the  center  of  gravity  of 
the  segment. 

The  procedure  for  the  experimental  determination  of 
the  mass  moment  of  inertia  matrix  using  either  a  torsional 
or  planar  pendulum  oscillation  technique  has  been  de¬ 
lineated  by  Wells  (4)  and  extensively  reviewed  and  de¬ 
veloped  by  Chandler  et  al-  (i). 

In  the  present  study  we  determined  the  center  of  gra  v¬ 
ity  (CG)  of  various  helmets  using  a  modification  of  the 
differential  weighing  technique.  The  torsional  pendulum 
technique  was  adapted  for  the  measurements  of  the  mass 
moment  of  inertia  of  helmets  about  the  CG  with  one 
plane  of  symmetry.  Having  determined  the  CG,  the  co¬ 
ordinates  of  any  point  on  the  rigid  body,  with  respect  to 
an  orthogonal  system  about  the  CG,  is  very  often  a  ne¬ 
cessity.  We  have  developed  an  electronic  caliper  for  length 
measurements  together  with  the  appropriate  software  to 
convert  these  length  data  into  three-dimensional  coor¬ 
dinates. 

MATERIALS  AND  METHODS 

Determination  of  the  inertial  and  geometrical  prop¬ 
erties  was  accomplished  in  three  stages;  1)  localization 
of  the  center  of  gravity  (CG)  and  placement  of  a  local 
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body  fued-Cartesian  coordinate  system  it  the  CG,  2)  de¬ 
termination  of  geometry  of  any  point  on  the  rigid  body, 
and  3)  determination  of  the  principal  mass  moments  of 
inertia  and  its  associated  directions  about  the  CC.  Each 
of  these  topics  will  be  addressed  individually. 

Center  of  gravity  determination.  A  modified  version 
of  the  differential  weighing  method  was  used  for  locali¬ 
zation  of  the  helmet’s  CG.  The  use  of  the  balance  plate 
for  differential  weighing  rests  theoretically  on  the  prin¬ 
ciples  of  statics.  The  experimental  arrangement  consisted 
of  an  aluminum  balance  plate  supported  by  knife  edges; 
the  supports  were  maintained  parallel  to  each  other  at  a 
known  distance  apart  by  placing  them  in  machined  grooves 
at  each  end.  Under  one  knife  edge  a  weight  balance  or 
a  force  transducer  was  mounted.  With  only  the  balance 
plate  in  place,  the  summation  of  moments  about  knife- 
edge  A  in  Figure  ia  yields: 

W*D*  -  FfD«,  Ill 

where  is  the  weight  of  the  balance  plate  acting 
through  its  CC  at  a  distance  Dt,  from  A.  and  F|  is  the 
initial  transducer  output  with  a  span  length  of  Dm  from 
A.  With  the  helmet  on  the  balance  plate,  we  get: 

WJ>,  +  W*D*  -  (F|  +  AF»)D»,  [2j 

where  Wk  is  the  helmet  weight  acting  through  its  CC  at 
a  distance  of  Dk  from  A,  and  AFa  is  the  increase  in  the 
transducer  output.  Subtracting  equation  1  from  2,  we  get 
either 

Dm  -  AFgWJD,  [3a] 

or 

AF,  -  WkIVD,.  13b} 

Because  WK  and  D*  can  be  predetermined  and  AF»  mea¬ 
sured.  one  can  either  mark  the  helmet  by  moving  some 
indicator  through  the  distance  Dfc  calculated  from  equa¬ 
tion  3a  or  by  fixing  Dk  and  moving  the  helmet  until  a 
calculated  AF«  from  equation  3b  is  read  on  the  trans¬ 
ducer. 

In  our  setup,  we  surrounded  the  balance  plate  on  three 
sides  with  a  C -channel  frame  positioned  such  that  it 
spanned  the  midsection  of  the  balance  plate  (see  Figure 
lb).  Across  the  top  of  this  frame  was  a  linear  series  of 
small  holes.  The  frame  was  placed  such  that  the  vertical 
plane  passing  through  these  holes  intersected  the  balance 
plate  along  a  line  parallel  to  the  knife  edges.  Care  was 
taken  to  ensure  that  the  vertical  plane  and  the  balance 
plate  surface  were  perpendicular.  Once  this  orientation 
had  been  achieved,  the  frame  and  plate  supports  were 
rigidly  fixed. 

A  Cartesian  coordinate  system  based  on  the  CG  was 
mapped  onto  the  helmets  using  the  following  procedure. 
Prior  to  mapping,  the  system  orientation  was  defined  as 
the  helmet  in  an  upright  position,  i.e.,  point  contact  with 
a  flat  surface  occurs  at  the  back  and  at  the  ear  flaps.  For 
those  helmets  where  these  points  were  poorly  defined, 
20-gauge  hypodermic  needles  were  placed  in  the  helmet 
shell  such  that  point  contact  occurred.  In  all  helmets,  the 
contact  points  on  the  ear  flaps  were  essentially  symmetric 
about  the  midsagittal  plane. 


the  V  of  the  helmet  ami  one  at  each  oar  flap. 


Figure  III  Iwuninir  rhw  of  C-ehannet  frame  surrounding  the  balance 
plate.  Sena  the  top  of  this  frame  are  a  linear  aerie*  of  small  boles  to 
aHow  a  pendulum  boh  carrying  a  felt-tip  pen  to  marli  the  outside  surface 
of  the  helmet. 

The  distance.  D*.  was  set  to  be  the  horizontal  distance 
from  edge  A  to  the  vertical  plane  passing  through  the 
holes  in  the  channel  frame.  Using  equation  lb.  a  value 
of  AFt  was  calculated.  The  helmet  was  moved  alone  the 
plate  in  planar  motion  until  the  calculated  AF(  coincided 
with  the  output  from  the  transducer.  At  this  point,  the 
CG  of  the  helmet  is  contained  in  the  vertical  plane  passing 
through  the  C-channel  frame.  The  helmet  was  then 
marked  using  a  pendulum  with  an  inked  felt-tip  pen  which 
was  successively  suspended  from  each  hole  in  the  frame, 
and  a  line  was  drawn  through  the  marked  points.  The 
helmet  was  then  routed  through  increments  of  approx¬ 
imately  60’  until  three  intersecting  lines  were  drawn.  A 
perspective  view  of  the  helmet  axis  system  and  inter- 
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secting  mapping  tines  is  shown  in  figure  2.  This  inter¬ 
section  ideally  should  be.  and  in  fact  was,  consist endy 
within  1  mm  of  being  a  point.  This  point.  Z,  on  the  helmet 
should  be  along  a  vertical  tine  passing  through  the  CG, 
Le.,  the  z-oxis  intercept. 

To  obtain  the  y-axis  intercept,  the  helmet  was  rotated 
until  the  z-axis  was  in  a  horizontal  plane.  This  was  ac¬ 
complished  with  the  aid  of  a  precision  rectangular  Plex¬ 
iglas  block.  The  block  was  placed  on  a  horizontal  surface 
and  the  helmet  positioned  with  the  aid  of  lightweight 
styrofoam  wedges  such  that  the  two  ear-flap  contact  points 
were  vertical,  i.e.,  along  the  edge  of  the  block.  Once  this 
position  was  attainted,  the  wedges  were  bonded  to  the 
helmet  surface.  The  differential  weighing  procedure  was 
repeated  to  obtain  the  + y-axis  intercept.  The  —y-axis 
intercept,  i.e.,  ~  Y,  was  mapped  with  the  aid  of  a  sharp 
pointer  inserted  into  the  center  of  the  balance  plate  at  a 
distance  of  D*  from  A.  The  helmet  was  placed  on  the 
balance  plate  with  the  +Y  point  on  the  pointer.  The  hel¬ 
met  was  then  tipped  with  the  aid  of  styrofoam  shims  until 
the  transducer  output  reached  the  calculated  AFt.  The 
plumb  bob  was  used  to  mark  —  Y  point.  The  x-axis  in¬ 
tercept.  X.  was  found  by  positioning  the  helmet  in  the 
reference  orientation  and  measuring  the  vertical  distance 
from  the  ear  flap  contact  points  to  ±Y.  This  measured 
distance  was  marked  as  a  line  above  the  contact  point  of 
the  hack  of  the  helmet.  The  x-z  plane  was  inscribed  by 
routing  the  helmet  forward  (about  a  line  parallel  to  the 
Y-axis)  about  the  ear-flap  contact  points.  The  intersection 
between  the  previously  marked  horizontal  line  and  the 
inscribed  line  determined  the  X  point.  The  distances  XY 
and  —  YX  should  be  equal. 


Upon  obtaining  the  three  marked  points  on  the  surface 
of  the  helmet,  the  actual  location  of  the  CG  is  found  by 
assuming  an  orthogonal  coordinate  system  which  inter¬ 
sects  the  three  points  simultaneously.  Figure  3a  illus¬ 
trates  the  geometrical  relationship  between  the  three 
points,  X,  Y,  and  Z,  and  the  CG.  By  the  Pythagorean 
Theorem,  we  get  three  nonlinear  algebraic  equations: 


a*  +  b*  -  A* 

141 

a*  +  c*  -  C* 

151 

b*  +  c*  -  8* 

[61 

where  A,  B,  and  C  represent  the  physical  distances  XY. 
YZ,  and  XZ,  respectively,  on  the  helmet  and  a,  b,  ami  c 
represent  the  distances  from  the  CG  to  the  points  X  Y, 
and  Z,  respectively.  Subtracting  [6]  from  [5}  and  adding 
the  resultant  to  |4j  yields 

a  -  1(C*  -  B*  +  A*)/2]*«.  17j 

Similarly, 

b  -  [{A*  ~  C*  +  B*)/2|1'*  18] 

c  -  [{8*  -  A*  +  C*)/21‘*  191 

where  positive  signs  were  used  in  the  above  equations 
due  to  physical  reasons. 


Figure  3a— The  gaootetncaj  relationship  between  the  canter  of  gravity 
(CG)  of  the  rigid  body  end  the  t,  y,  and  z  axis  intercepts!  X.  V,  and  Z. 
at  well  as  a  fanatic  point  (f ). 
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Figure  3b— The  pmfttetritel  relationship  between  the  center  of  gravity 
(CG)  of  tt>e  tipd  body  and  the  s  V  intercepts  as  weii  as  a  generic  point 
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Tbree-dimeauoujtl  coordinate  determination.  Having 
an  orthogonal  coordinate  system  fixed  to  the  CG,  it  is 
theoretically  possible  to  obtain  the  coordinates  of  any 
generic  point,  P,  on  the  helmet  Referring  again  to  Figure 
3a,  we  note  from  the  law  of  cosines  that 

L*  -  K*  4-  a*  -  2aR,  [10a] 

M*  -  R*  +  b*  -  2bR,  [10b] 

N*  -  R*  +  e*  —  2cR*.  [10c] 

where  R„  R,,  and  R,  represent  the  vector  components  of 
>S),  the  position  vector  from  the  CG  to  the  point,  P.  Be¬ 
cause  the  direction  cosines,  l  ns,  and  n,  of  the  vector  R 
are  defined  by  :  1  —  RJfRi,  m  *•  R,®!,  and  n  *  RJlRl, 
we  obtain  R,.  R,,  and  R,  from  equation  10  to  yield 

I  -  (a*  -  R*  -  L*)/2!R_ia  [llaj 

m  -  (b*  +  R*  -  M*)/2iRlb  (llbi 

n  -  (c*  +  R*  -  N*)/2&te.  [lie] 

In  theory,  if  we  can  measure  the  distances,  L,  M,  N,  and 
®1.  we  can  solve  for  die  direction  cosines  via  equation 
1 1  because  a,  b,  and  c  are  given  in  [7]  —  (9|.  Practically, 
however,  R1  is  physically  difficult  to  obtain.  To  get  iRl. 
we  take  advantage  of  midsagttta!  plane  symmetry  and 
pick  a  convenient  auxiliary  point,  —  Y,  on  the  shell.  Re¬ 
ferring  to  Figure  3b,  we  apply  the  law  of  cosines  twice: 
once  for  the  entire  triangle  and  Oi.ce  for  the  smaller  tri¬ 
angle  on  the  right,  it  can  be  easily  shown  that 

R*  -  (M*  +  K*  ~2b*)/2,  [12] 

where  K  is  the  distance  from  — Y  to  the  generic  point. 
P.  Thus,  by  making  an  additional  distance  measurement, 
we  can  obtain  R. 

In  practice,  we  used  equation  12  to  estimate  R,  sub¬ 
stituted  it  into  equation  1 1  to  obtain  the  direction  cosines 
l,  m.  and  n.  These  calculated  values  were  substituted  into 
the  well-known  constraint  equation: 

1*  +  m*  +  a1  -  I.  [13] 

If  the  sum  of  the  squares  of  the  direction  cosines  t$  l  ± 

« (where  t  is  a  small  number),  the  values  of  i,  m,  and  n 
were  considered  adequate.  Otherwise,  the  value  of  R  was 
iterated  until  the  criterion  was  met.  The  above  iteration 
procedure  is  in  the  form  of  an  efficient  computer  algo¬ 
rithm  and  yields  results  almost  as  soon  as  the  physical 
measurements  are  completed.  Thus,  given  the  distances 
A,  B,  C,  L,  M.  N,  and  K,  a  simple  vector  analysis  yields 
the  position  vector,  i.e.,  three-dimensional  digitization, 
from  a  reference  point.  What  is  required  is  a  device  to 
accurately  measure  these  distances. 

Electronic  caliper.  An  electronic  caliper,  shown  in  Fig¬ 
ure  4.  was  designed  to  make  these  linear  measurements. 
It  consists  of  a  three-bar  planar  linkage  arranged  in  a  U- 
shape  configuration.  The  two  hinge  joints  are  two  linear, 
high-precision  rotational  potentiometers  The  voltage 
output  of  each  potentiometer  is  proportional  to  the  angle 
between  the  associated  free  arm  and  a  known  reference 
position.  The  analog  signals  are  converted  to  digital  ones 
through  a  minicomputer  (PDP  11/34)  and  all  subsequent 
processing  is  in  the  digital  mode.  Digitization  is  at  a  rate 
of  1  kHz  and  duration  of  0.1  s,  A  foot  switch  triggers 
A/D  conversion  and  enters  the  distance  being  measured 
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Figure  4— Tto  artaccraaic  c*Uper  urod  to  obtain  length  mtattuanauls. 
Tin  two  hinge  joint!  ant  two  linear,  high-precision  rotational  potto- 
tlonttaa  The  voitagt  output  of  each  potentiometer  it  proportional  to 
that  ingie  Iwtwwtn  the  free  arm  tad  a  known  reference  paaitkm. 


by  the  hands  holding  onto  the  two  arms  of  the  caliper. 
Accuracy  was  ±  0.4  mm  over  a  range  of  40  cm. 

Using  the  caliper  device  an  operator  can  do  three- 
dimensional  digitization.  In  this  mode  the  first  three  dis¬ 
tances  entered  determine  an  orthogonal  coordinate  sys¬ 
tem.  Each  subsequent  set  of  four  distances  uniquely  de¬ 
termines  the  Cartesian  coordinates  of  a  given  point  relative 
to  the  system  established.  As  stated  previously,  the  value 
of  R  was  iterated  because  its  input  value  was  only  ap¬ 
proximate.  This  iteration  was  continued  until  the  right 
side  of  equation  13  was  equal  to  1  £  0.01.  Convergence 
was  always  achieved  within  30  steps  and  usually  within 
10.  While  the  electronic  caliper  was  used  to  determine 
the  coordinates  of  selected  points  on  the  helmet  shell  in 
this  project,  it  is  potentially  suitable  for  three-dimen¬ 
sional  digitization  of  any  geometry. 

Mass  moment  of  inertia.  If  one  assumes  midsagittal 
plane  symmetry,  e  g.,  x-z  plane,  the  mass  moment  of  in¬ 
ertia  matrix  about  the  CG  reduces  to  five  non-zero  com¬ 


ponents: 


I„  0  1., 

li  =“  <>  L  <> 


The  method  for  the  determination  of  the  mass  moments 
of  inertia  about  a  given  axis  via  the  torsional  pendulum 
technique  is  well-known,  see  e.g..  Liu  et  ai.  (2).  The  oniv 
non-zero  off-diagonal  term  is  I„.  It  is  quite  well-estab¬ 
lished.  see  e.g..  Wells  (4).  that 
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where  is  the  mass  moment  of  inertia  about  an  arbitrary 
q-q  axis  and  X.  m.  *  are  its  direction  cosines. 
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Flgura  ft-Tb*  nan  of  Army  aviator,  Aaaariean  football  ami  bicycle 
Wmett-  The  abbveviatwoa  rap.wurat  tb*  followin®  I)  Bnutesie— bicycle 
beiaM*;  31  Bid  Kralitv- BiddeS  football  helmri  Model  Kndite  3)  Army 
A-l-Army  aviator  exp— ia—otai  helmet;  4}  Bifcc  IV-B&e  football  W- 
aaat  Modal  IV;  5)  Raw.  HC— RawUagi  football  helmet  Model  HC;  6) 
Sid.  PAC^d-  Riddell  football  helmet  Modal  T AC-44-,  7)  Mac. 
106MH»  MacGregor  football  helmet  Model  100-MH;  8)  Raw. 
HKEhP-Hawiiam  football  halowt  Model  KND-P;  9)  Bid.  ACE-1  -RWklell 
football  bab»«t  Modal  ACE-l;  18)  Wtt.  F-tlOl  -  WiLoa  football  helmet 
Medal  FtlOJ;  U)  Army  SPH-t— Army  aviator  bebaet  Modal  SPH-t; 
end,  IS)  Saw.  KNOd-Sawbaga  football  helmet  Model  HND-9. 


Experimentally,  was  obtained  via  the  torsional  pen¬ 
dulum  technique  and  the  direction  cosines  A,  m.  and  > 
were  found  by  using  the  electronic  caliper  and  the  com¬ 
puter  algorithm. 

To  obtain  the  principal  mass  moments  of  inertia  and 
its  associated  directions,  we  can  use  either  Mohr's  circle 
or  classic  analysis  (3).  In  either  case,  one  gets 
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where  It,  It,  Ij  are  the  principal  mass  moments  of  inertial 
and  8  is  the  needed  rotation  of  the  x  and  z  axes  about  the 


y-y  axis  to  yield  the  principal  directions  as  shown  in  Fig¬ 
ure  5. 


of  the  helmet,  as  the  reference  orientation,  the  CG  is 
located  by  a  vertical  distance.  Zc.  from  a  flat  surface.  The 
calculated  distances  a.  b,  and  c  serve  as  a  measure  of 
helmet  size.  These  results  are  given  in  Table  l.  From  this 
reference  orientation,  the  distance  from  the  CG  to  any 
point  can  be  measured.  For  example,  when  worn,  the 
distance  from  the  helmet's  CG  to  the  zero  tangent  point 
of  the  ear  flap  can  be  easily  determined  using  the  3D 
digitization  scheme  delineated  earlier. 

The  product  moment  of  inertia,  i.e..  I„,  and  the  angular 
rotation,  8,  to  obtain  the  principal  directions  are  shown 
in  bar  graphs  in  Figure  7.  The  mass  moments  of  inertia 
about  the  orthogonal  coordinate  system  placed  at  the  CG 
of  the  helmets  in  the  reference  configuration,  with  and 
without  a  face-guard,  are  shown  in  Table  2.  The  principal 
mass  moments  of  inertia  together  with  the  angle  of  ro¬ 
tation,  9,  defining  the  principal  directions  are  given  in 
Table  3. 


RESULTS 

The  weight  (mass),  center  of  gravity  (mass),  and  the 
principal  mass  moments  of  inertia  of  one  bicycle,  two 
Army  aviator,  and  nine  American  football  helmets  were 
obtained.  The  mass  in  grams,  with  and  without  the  face- 
guard,  are  shown  in  Figure  6.  Using  the  upright  position 


DISCUSSION 

Prior  to  the  CG  and  the  mass  moment  of  inertia  de¬ 
termination,  the  instruments  were  all  calibrated.  The  CG 
determination  was  checked  with  a  high-precision  ma¬ 
chined  Plexiglas  rectangular  block.  Its  x-,  y-,  and  z-axss 
intercepts  are  easily  determined  by  geometry  and  checked 
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Figure  7— The  product  moaeut  of  inertU,  t^,  erd  the  uguler  rotation. 
#,  retjuired  to  obtain  the  principal  directions,  usurping  otidstgittai  plane 
jymrnetry,  for  the  hdtmti  studied. 


against  the  balance  plate  method.  The  resolution  of  the 
CG  determination  technique  is  0.5  mm.  The  location  oi 
the  x-  and  z-axis  intercepts,  i.e.,  X  and  2  points,  invariably 
coincided  closely  with  mold  lines  for  the  helmet  shell, 
affirming  the  essential  correctness  of  ‘he  midsagittal  plane 
symmetry  assumption. 

The  3D  digitization  method  was  compared  with  one 
that  has  been  commercially  available  (Oral-Pen  3-D  Dig¬ 
itizer,  Science  Accessories  Corporation,  Southport,  CT 
05490).  The  operation  of  the  Graf-pen  system  utilizes 
the  essential  constancy  of  the  speed  of  sound  in  air.  A 
spark  gap,  placed  at  the  point  ;o  be  digitized,  when  Bred, 
creates  a  spherical  sound  wave  which  is  picked  up  by 
thr-i-e  linear  microphones  arranged  as  an  orthogonal  ar¬ 
ray.  Ti  e  time  it  takes  the  sound  wave  to  reach  a  given 
microphone  yields  the  magnitude  of  the  distance  given 


the  constancy  of  the  velocity  of  sound.  The  system  has 
the  advantage  of  speed  but  is  not  without  its  problems, 
the  spark  gap  sometimes  misfires,  and  there  are  restric¬ 
tions  in  terms  of  the  location  of  the  spark  gap.  No  ob¬ 
structions  are  allowed  between  the  spark  gap  and  the 
linear  microphones.  When  compared  against  this  sophis¬ 
ticated  and  expensive  system  on  a  standard  rectangular 
calibration  block,  our  system  was  more  accurate,  even 
though  it  took  four  times  as  long.  On  experimental  hel¬ 
mets,  the  two  methods  provide  similar  results. 

The  center  of  gravity  (mass)  and  the  principal  mass 
moments  of  inertia  of  a  rigid  body  are  measures  of  the 
mean  and  the  distribution  of  mass  points  about  the  mean, 
i.e.,  about  an  orthogonal  coordinate  system  fixed  at  the 
center  of  gravity.  Front  the  theoretical  point  of  view,  the 
helmet  may  b  *  idealized  as  a  hollow  ellipsoid  whose  prin¬ 
cipal  mass  moments  of  inertia  about  its  CG  are: 

L  “  [m,(b*  +  e3)  -  m,(b?  4-  c?)]/5 

I„  -  fm.{a*  +  c3)  -  m,(af  4-  c?)]/5  (18| 

I„  “  (m„{a3  +  b3)  -  m,(a,2  +  b?)]/5, 

where  m0  and  m,  are  the  mass  of  outer  and  inner  ellipsoid, 

respectively,  and  a,,  b„  and  c,  are  the  inner  ellipsoid  axis- 
intercepts  with  the  x,  y,  and  z  axes,  respectively. 

If  one  defines  the  characteristic  thickness  in  the  three 
coordinate  directions  as: 

r  -■  a  —  a, 

s  “  b  ~  b,  [  19| 

t  *  c  -  c, , 

it  is  easily  shown  that 

_  L  =  lm(b*  +  c*)/5{  +  m,[(b*  +  c’)  -  (b?  4-  c?)}/5 
4,  **  [m(a3  4-  c!)/5j  4-  m,!(aJ  4-  c1)  —  (a?  4-  [20| 

I„  =»  [m(a*  4-  b*)/5|  4-  m,l(a3  4-  b3)  -  (a?  +  bf)]/5. 
where  in  ”  m„  —  m,  ~  the  mass  of  the  hollow  ellipsoid. 
Thus,  for  the  nonhomogeneous  helmets  considered  here, 
the  most  appropriate  form  of  the  regression  equation  for 
the  three  principal  mass  moments  of  inertia  are: 
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I«  -  Cim(b*  +  c»)  +  C*(b*  4-  c1)  +  C, 

I„  -  C,m(a‘  +  c*)  +  C,(a*  +  c*)  +  C,  (211 

L  -  0,111(4*  +  b»)  +  C,(a*  +  b*)  +  C* 
where  C,,  C*  .  .  .  C,  are  regression  coefficients. 

For  the  helmets  studied,  it  appears  that  the  weight 
(mass)  and  the  geometrical  distances,  a,  b,  and  c,  are  highly 
correlated  with  the  principal  mass  moments  of  inertia  as 
shown  in  Table  4.  For  instance,  even  when  football  and 
Army  aviator  helmets  were  mixed,  the  correlation  coef¬ 
ficients  for  the  linear  regression  ranged  from  0.88  to  0,98 
in  Table  4. 

To  validate  these  regression  equations,  a  Riddell  Model 
PAC-3  football  helmet  was  weighed.  Substituting  its  mass 


of  1049.9  g  into  the  regression  equations  for  football  hel¬ 
mets  without  face-guards  in  Table  4  yielded  the  predicted 
inertial  properties  shown  in  Table  5.  The  difference  be¬ 
tween  experimentally-determined  and  the  predicted,  val¬ 
ues  shown  in  Table  5  is  less  than  3%. 

An  examination  of  Table  1  showed  that  a,  b,  and  c  for 
football  helmets  are  all  in  the  range  of  10.5  to  12.6  cm. 
i.e.,  an  ellipsoid  with  eccentricities  small  enough  to  be 
considered  practically  a  sphere.  For  a  sphere  of  radius 
R,  the  mass  moment  of  inertia  about  any  diameter  is  I  “ 
2mR*/5.  It  is  not  too  for-fetched  to  assume  that  the  prin¬ 
cipal  mass  moments  of  inertia  might  be  linearly  correlated 
with  mass,  i.e.. 


TABLE  2.  Tht  numerical  values  of  the  mm 
and  without  s  lace-guard.  Midsagittal  plane 


martu  ahout  ttM  erthogeoM  eoordiwrtt  system  tad  to  the  CG  ot  the  varans  helmets  in  the  referent*  configuration  with 
of  the  helmet  it  assumed. 
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TABLE  3.  The  numerical  value  of  the  principal  mass  moments  of  inertia  together  with  ttM  angular  rotation,  S,  required  ^  y*id  the  principal  directions.  Midsagittal  plane  symmetry 
is  tssumed. 
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dtotrmmd  Mat 


CifcotaM 

rmaf* 

tipiriawatif 

rrwrtpw 

* - al-- 

fssnxn 

PwtseSa** 

Error 

ftddtU 

1*  -  133 

It  -  6.1 

2.68 

PAC-3 

iy  -  9.3? 

iy-  12 

1.92 

*  -  1049.9 1 

It  -  10  .31 

k  -  10.2 

1.03 

TABLE  6.  Tba  linear  rtgrtssioa  agaatwa  taUting  tot  priaeipai  aast  aoawits  to  iaartia 
to  tot  sun  to  football  ktoaati 

Uaaar  Bagrassiaa 

Camtattta 

Cf  rfftiM 

*ftadpai  laartiat 

CaafBdaat 

Casntanto  Amy  and 

lt  -  0.0078(a)  +  0.459$ 

0.7502 

Football  Holatatt 

I,  -  0.0102<m)  -  1.6414 

9.8011 

(w/c  fuanli) 

l?  -  0.0107(a)  -  0.7238 

0J866 

1,  -  0.0077(a)  +  0.2366 

0.9261 

Football  Hstmtos 

1,  -  0.0119(a)  -  3.0531 

0.9768 

(w/o  punts) 

4  -  0.0132(a)  -  12626 

0.9332 

I  “  Cl#in  +  C„,  [22] 

where  C,»  and  C,  t  are  linear  regression  coefficients.  Table 
6  yields  the  linear  regression  equations  for  the  three  prin¬ 
cipal  mass  moments  of  inertia  for  the  helmets  studied. 
Note  that  the  correlation  coefficients  are  only  slightly 
lower  than  those  in  Table  4.  Figure  8  displays  the  ex¬ 
perimentally-determined  principal  mass  moments  of  in¬ 
ertia  as  a  function  of  the  football  helmet  mass  together 
with  their  linear  regression  lines.  The  slopes  of  these 
regression  lines  are  slight,  i.e.,  for  the  football  helmets 
examined  the  inertial  properties  increase  slowly  with  mass. 
For  instance,  the  slope  of  i,  in  Figure  8  is  0.013  g-mJ/g. 

If  one  were  to  substitute  the  Riddell  Model  PAC-3  data 
into  the  linear  regression  equations  shown  in  Table  6,  the 
difference  between  the  experimentally-obtained  data  and 
the  predicted  values  is  less  than  4%.  Thus,  in  any  dynamic 
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Figaro  8 — The  linear  regression  equation,  relating  the  principal  mats 
moments  of  inertia  to  the  mass  of  the  helmet  together  with  the  exper¬ 
imentally-determined  principal  mass  moments  of  inertia  of  actual  foot¬ 
ball  hairnets  without  face-guards. 


Study  involving  the  current  generation  of  football  and 
Army  aviator  helmets,  one  can,  given  their  mass,  estimate 
their  principal  mass  moments  of  inertia  to  within  ap¬ 
proximately  5%  of  their  experimentally-determined  val¬ 
ues.  The  standard  error  of  estimate  for  principal  inertias 
of  all  helmets  used  in  this  study  was  0.538  g-mJ  when 
using  the  linear  regression  equation.  The  above  results 
are  possible  because  the  mass  distribution  of  the  current 
generation  of  helmets  is  primarily  set  by  head  size  and 
secondarily  by  helmet  size,  shape,  and  materials.  Should 
the  size,  shape,  and  material  arrangement  in  the  next 
generation  of  helmets  be  altered  substantially,  the  more 
complex  regression  equations  given  in  Table  4  will  be 
more  appropriate. 

The  praaant  project  was  principally  supported  by  the  Htoan  Streiftar 
Fund  of  tha  University  of  Iowa  Foundation.  Partial  funding  was  provided 
by  Contraa  No.  DAMO  17-8t-C-1 186  from  the  U  S.  Army  Aeromedica! 
Research  Laboratory.  Ft.  Rucker,  Alabama  and  Grant  No.  GM  26608- 
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Appendix  B 


The  node  shapes  of  SPH-4  (Fig.  Bl>  and  A1  (Ftg.  B2)  Army  Aviator  Helmets, 
eigenvalues  corresponding  to  these  mode  shapes  are  given  in  Table  2. 
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